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Целью работы является построение и исследование адаптивного 
метода минимальных поправок (МАММП), предназначенного 
для численного моделирования процессов фильтрации двухфаз- 
ной сжимаемой жидкости в пористых средах. Описанный под- 
ход позволяет преодолеть известные ограничения использова- 
ния других методов решения сеточных уравнений, определяе- 
мые такими условиями, как: значительные перепады давления, 
действующего на нефтеводоносный пласт; сжимаемость среды 
при значительном содержании газа в нефтяной фазе.В качестве 
основы исследования выбран метод аппроксимации — явный 
для определения функции водонасыщенности и неявный для 
расчета функции давления.При постановке начально-краевой 
задачи и ее дискретизации рассмотрен процесс фильтрации 
сжимаемой  двухфазной жидкости в 
трехмерной области с боковой поверхностью, ограниченной 


пространственно- 


снизу поверхностью подошвы пласта, а сверху — поверхностью 
его кровли. Построен двухслойный итерационный метод вариа- 
ционного типа — модифицированный метод минимальных по- 
правок, адаптированный для решения сеточных уравнений 
двухфазной сжимаемой жидкости с несамосопряженным опера- 
тором при самых общих предположениях относительно свойств 
оператора сеточной задачи. Показано, что МАММП обладает 
асимптотической скоростью сходимости, характерной для 
«классического» ПТМ, не использующего технику чебышевско- 
го ускорения и применяемого для задач с самосопряженным 
оператором. Численные эксперименты подтвердили высокую 
эффективность МАММП. Установлено, что для достижения 
заданной точности число итераций при МАММП сокращается в 
3-20 раз по сравнению с методом Зейделя и методом верхней 
релаксации. 


Ключевые слова: задачи фильтрации двухфазных сжимаемых 
жидкостей, сеточные уравнения, несамосопряженный оператор, 
адаптивный метод минимальных поправок. 
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Введение. Адаптивный метод минимальных поправок (МАММП) предназначен для численного моделирова- 
ния процессов фильтрации двухфазной сжимаемой жидкости в пористых средах. Такие процессы являются базовыми 
в задачах проектирования разработок нефтяных месторождений в ряде случаев. Во-первых, если нефтеводоносный 
пласт находится под воздействием давления, изменяющегося значительным образом (в несколько раз или даже на 
один-два порядка). Во-вторых, если нефтяная фаза имеет значительное содержание газа и необходимо учитывать 
сжимаемость среды. В этих условиях дифференциальный оператор уравнения для определения давления становится 
несамосопряженным, что приводит к несамосопряженности оператора сеточной задачи. Для задач фильтрации сжима- 
емой жидкости в самой общей постановке также не является ограниченным сеточное число Пекле. Эти особенности 
задач фильтрации двухфазной сжимаемой жидкости существенно ограничивают возможности использования эффек- 
тивных методов решения сеточных уравнений (таких, как предложенный авторами ранее усовершенствованный попе- 
ременно-треугольный метод, ПТМ). В качестве примеров можно привести численное решение задач фильтрации не- 
сжимаемой жидкости с самосопряженным оператором, а также уравнений с несамосопряженными операторами и 
ограниченным значением сеточного числа Пекле. 

В данной работе, как и в [1], авторы ориентируются на метод аппроксимации — явный для определения 
функции водонасыщенности и неявный для расчета функции давления (1МРЕ5-метод) [2]. В этом случае в построении 
и численной реализации дискретной модели для определения функции давления возникают дополнительные трудно- 
сти, связанные с изменением плотности более сжимаемой нефтяной фазы. Это, в свою очередь, приводит к изменению 
свойств оператора уравнения — он становится несамосопряженным [3], [4], и в общем случае невозможно гарантиро- 
вать ограниченность сверху величины сеточного числа Пекле, например, константой 2. Поэтому ни один из вариантов 
ПТМ, построенного в работах [1], [5], [6], не пригоден для численной реализации модели. Сравнение построенного 
метода с употребительными подходами к решению сеточных уравнений для восстановления функции давления (метод 
Зейделя, метод верхней релаксации) показывает существенное преимущество модифицированного адаптивного мето- 
да минимальных поправок (МАММП). Выигрыш в числе итераций составляет 3—40 раз в зависимости от числа ячеек 
используемых сеток, что позволяет использовать МАММП в качестве базового алгоритма решения сеточных задач 
проектирования разработок нефтяных месторождений. 

1. Постановка начально-краевой задачи и ее дискретизация. Как ив [1], рассмотрим процесс фильтрации 
сжимаемой двухфазной жидкости в пространственно-трехмерной области С с боковой поверхностью 


5(х,у), (ху)еГ, Г=В\ О, ограниченной снизу поверхностью подошвы пласта 


&=Нь (х, У), (х, У) ЕО, ОР ЕК”, 
а сверху — поверхностью кровли пласта 
<=Н, (х, У), (х, у) ЕЛ. 
Запишем уравнения неразрывности для каждой из двух фаз в отдельности в дивергентной форме. 
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где индексы = 1,2 обозначают соответственно нефтяную и водную фазы; р;, 5; — плотность и насыщенность соот- 


р 
— 
ветствующих фаз; и; — скорость их фильтрации; 4; — мощность объемных источников (стоков) каждой из фаз. 


Будем ориентироваться на закон Дарси, который запишем в виде уравнений для каждой из фаз в отдельности: 
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где К — коэффициент абсолютной проницаемости; }; = 2 (52) — относительная проницаемость соответствующей 
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Уравнения (1) и (2) дополняются уравнениями состояния для каждой из фаз: 


р; = ры [1+В:(р-2)], #=12, (3) 


т 
> > > д д д >> 
фазы; в — вектор ускорения свободного падения; Квт = [№ Е РК, 7. ‚ К. 8 = (0,0, 8)". 
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где ро;› Ро — соответственно характерные значения плотности и давления, В; — коэффициенты сжимаемости сред. 


Предположим, для определения функций }1(52) и №(52) используются многочлены третьего порядка относи- 
тельно функции водонасыщенности 52 [2]. Чтобы получить уравнение для определения давления, подставим в каждое 
из двух уравнений (1) выражение для скорости фильтрации соответствующей фазы из уравнений (2). Таким образом, 
придем к соотношениям 
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Разделим уравнение (4) на произведение функций тр, а уравнение (5) — на тр» с учетом уравнений (2). 
Складывая преобразованные уравнения, получим: 
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При этом условии соотношение (7) является уравнением параболического типа, содержащим в общем случае 
младшие производные (первого порядка по пространственным переменным) относительно плотности. Учитывая (3), 
уравнение содержит первые производные относительно давления, а также переменные множители вида 
Итр;, 1=1,2 при операторах 

> > > 
@» ВЫ бро :] ‚ 1=12. 
Ш 

Итак, дифференциальный оператор в правой части уравнения (6) является несамосопряженным [3]. Получена 
система выражений, включающая уравнение (5) или вместо него — (4), уравнение (6), уравнение состояния (3). К 
данной системе необходимо добавить соответствующие начальные и граничные условия, сформулированные в [1], а 
также уравнение (4). Для численной реализации данной модели целесообразно провести линеаризацию задачи (по 
временной переменной), выделить в операторе задачи на верхнем временном слое симметричную и кососимметрич- 
ную части, а затем выполнить дискретизацию. Построим неравномерную временную сетку: 


6. = в ь тк, 15изм} 
К=1 


Для водонасыщенности $2 опустим нижний индекс и построим явное разностное уравнение (по времени). Бу- 
дем считать пористость т не зависящей от давления. 
Из уравнения состояния (3) получено равенство 


5: Ар РИ Чр› _ (1-5) В, р: 5В> | 
р 4 р› 4 1+В. (р-р) 1+В> (Р-Ро) 





(8) 


Преобразуем первое слагаемое в правой части уравнения (7), выполняя линеаризацию на временной сетке @, . 
При этом отнесем линеаризованные коэффициенты на нижний временной слой (1=1, ), а члены, определяющие дав- 
ление, — на верхний временной слой. Для определенности {1 =1,.| в следующих ниже выражениях эти члены отмече- 
ны символом «/^» над функциями). 
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Аналогично предыдущему приходим к равенству: 


та о [45 ы | - 
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Подставляя (8), (9) и (10) в соотношение (6), получаем линеаризованное дифференциально-разностное урав- 
нение относительно давления с младшими производными: 
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Для функции водонасыщенности построим на временной сетке явное относительно 5 уравнение: 


т(р>$-—р25) а ие [4 о =] т 


Ти | т 





которое можно переписать в развернутом виде: 
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Для вычисления плотностей водной среды на нижнем и верхнем временных слоях используем уравнения со- 


стояния вида (3): 
Б> = Ро> [1+В» (Р- в») |, 
р> = Ри: [1+Во (р- рэ} |. 


Начальные условия для давления, которые требуются для решения уравнений (11) и (12) на первом времен- 
ном слое, рассчитываются, исходя из гидростатического приближения и известного (раздельного) распределения двух 
фаз до начала процесса фильтрации, с заданной границей раздела двух фаз. Кратко обсудим дискретизацию по про- 
странству полученных дифференциально-разностных уравнений (11), (12). Дифференциальные операторы, входящие 
в выражения вида: 


(13) 





| (-)в 5В> } 








+ + 
1+8 (Р-Р) 1+В>(р-ро))т 

+ ВРов 0 | Л(5) РИ Ро ВЕ 8 (5) @ , 4) 

мот 0” | ы тр; и; 0 





ый. | 1 (5) $ 2]. 4, | 
т| дх м; 0х) ду ш; ду} 05 и; 0 рут 


являются самосопряженными. Данное свойство сохраняется при использовании дискретизаций этих выражений на 
пространственных сетках — например, так, как это описано в [7]. Аппроксимация выражения (14) на сетке и одновре- 
менно определение самосопряженной части оператора сеточной задачи на сетке имеет вид: 

1 








р (1) р р (1) 
(^Р) ›=—8® (В д-Рд)В®. 
К й: тео 1] 
2 2 
ВФ (р. (>В а) в®, (Р.Р 
а и И РЕ ВА ыы] 
й; > й. 1,]--К 
и. 
2 Е и 
2 
о в 
АРА ыы] 1 РА РЬ] Вт 
й: А+ й: ВрЕ— 
(1-5 ле В 5; кВ т, к > 
+ + В+ 
я БЛ 
1+8, (Р БК - Ро) 1+8» (Р БК - Ро) п 
2 ВаРоа& ЕЕ = К Ле „к- | ^ 
+У В т. 
21. Иа 21. 
Х = я 
+к, Ро’ Ва 8 8 Зак В, ьы- Воды А й 
Ра: ‚.К На 21. 
0, 
+ Во, ы м. ‚рибь лк» если узел находится на нагнетательной скважине; (15) 
Ви, ‚Р ‚ к› если узел находится на эксплуатационной скважине. 
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При произвольных значениях шага по времени т. , а также при неизвестных заранее знаках производных 


2) „ 


д 
вида К, 
0 | 0% 


мирующего соответствующий дифференциальный оператор. Достаточное условие положительности сеточного опера- 
тора, аппроксимирующего дифференциальный оператор в соотношении (14), принимает вид: 


(1-5 ль В, + 5; кВ 


ое а 1+8, (рез — Ро) 1+В> (Рек — Ро) ое. 6) 


(коза 2 | Верь р дива] 





нельзя в общем случае гарантировать положительность сеточного оператора (15), аппрокси- 








а=| Й.т, 


‚„кНа Рав к 


Отдельно рассмотрим аппроксимации выражений в уравнении (11), содержащих производные первого поряд- 
ка относительно давления на верхнем временном слое: 


т Вгро: Л (5) Эр ФР | Вгро: Л (5) др. ФР, 
в й 
тр; |; Ох вх тр; |; 0» ду 











р (17) 
Ви Л) 12] > 
т | (р: / 0 
Обозначим 
р (ЕЕ, ) = ВаРоа о (5) др. (х,у, =) ЕС, 


Ра Ш 9 
В (х+0,5й, , У, СВ) = Вл (х +0,5й„Уу»льй, }, (ж‚Урль)е ог, 


ар(х,+0,51,,Уу»тььт, ) и. 








1 в 
® =©@ хо, хо, = - —, 18 
п 1 Хо Хоз 5х г (18) 
1 (х-0,5й,, У, 2,1, )= ь (х = УнЕь, ь. (журна) ор : 
Ор[х, —0,51., У; Зк> И, Р‚,-Р_ 
он = ХФ, Хх, ( В а В 1,2, 
Ох й 


х 





Е ВоРоа Ла (5) бр 
Й 
Рю Ш 9 


Ь, > (х,У+0,51., 2,1, = В. (У, +0,5й› 2», }, (х,у») =, 


БЕ = › (х,у, <) ЕС, 


Бык — Вл 
) 
й 


У 


2. _ о В 
®, =, Х0> Ход, = 





(19) 


В. 2 (х,У-0,5,,<.1,)= ь ( (ху, а (жур) ь 


@( ^»У)- г _ВоЕ-В 
= о | 


у 


2 
ре + 
9, = хо Хо, 





=1,2. 
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№(5)[ тд 
В (же унлие, = Е Вориа ( | (х,у,<) ЕС, 
а (Ра 63 


В. (х,у,=+0,51, 1, )= Ва (Уна +0, и, (%,у»)е вк, 


др(х,, У, <к +0,5й, „1 в... СР. 
9 = хо, хз, (х ] — ь) = ми м. (0) 





< 
В. (х,у,=-0,5,,1,)= Ы (жук — О |, (х,у» евр, 


др(х,, у, —0,5й., <. Р.к-Р; 
3 , ) 25К> ЗЕ 
[ум = вхо. хо, — 5% - ‚) ЕЕ. * ыы и а=12. 


д 
< 





С учетом введенных обозначений (18)—(20) сеточный аналог непрерывного выражения (17), умноженный на 


функцию т, ; ь принимает вид: 





а=1 


(АР) = К Вароа Л (5 ) бр Е, ВоРоа а $) бр бр, 
Ор р шщ бб рр ош 9 8 





Ра 0 05 | ом 


[и 


+ ВаРоа та, = Е (х+0,5й,,у,зи,)Р, + 
+61 (х-— 0,58, уг, )В;)+(6,, 2 (%у+0,5в, „зн, В, +62 (х,у-0,51, „2,1, ) В, + (21) 
+ (в. (х,у,=+0,5й.,1,)Р.+Ь, з (х,у,=-—0,5й, „1, ) В: ь.)], 


(мур зк) Е, хо, хоз. 


Таким образом, в операторном виде сеточную задачу, аппроксимирующую уравнение (11), можно записать в 


виде: 
АР=АР+АР= Х(Р,5), (22) 
А=А+А, А=А’, А=- 


Правые части определяются очевидным образом и отличны от нуля в узлах сетки, совпадающих со скважина- 
ми, а также в приграничных узлах сетки. Далее для простоты будем предполагать задание на границе области гранич- 
ных условий первого рода. 

2. Построение адаптивного метода минимальных поправок. Далее построен двухслойный итерационный 
метод вариационного типа — модифицированный метод минимальных поправок, адаптированный для решения се- 
точных уравнений двухфазной сжимаемой жидкости с несамосопряженным оператором [7] при самых общих предпо- 
ложениях относительно свойств оператора сеточной задачи. Требуется лишь, чтобы в явном виде были выделены 
симметричная и кососимметричная часть оператора сеточной задачи. Сформулируем задачу в общем виде. В конеч- 


номерном Гильбертовом пространстве Н рассматривается задача на обнаружение решения операторного уравнения: 
АБ=Г, А:Н->Н, (23) 


где А — линейный, положительно определенный оператор, определенный выше соотношениями (15), (21), (22),(А>0). 


Для нахождения решения задачи (23) используем неявный итерационный процесс: 


т-1 рт 


В +АР" =Е, В:Н>Н. (24) 
бт 





Здесь т — номер итерации; а. > 0 — итерационный параметр; В — легко обратимый оператор. Последнее 
означает, что обращение оператора В в (24) должно быть существенно проще, чем обращение исходного оператора А. 


В нестационарном итерационном методе вариационного типа 


рт-1 рт 
ВР ГАР" У (25) 


4-1 





итерационные параметры 01-1 выбираются из соображений минимизации скалярного произведения для погрешности 


решения в энергетическом пространстве Нрь [7]: 


Сухинов А. И. и др. Математическое моделирование фильтрации двухфазной сжимаемой жидкости 


^ 


5" = Р" _Р, $" = р"и _Р . (26) 
(Ре, =") 5 шш, р=Р’>о0. (27) 
Учитывая условие (26), а также представление ТН ато В-ТА2т ‚ приходим к задаче 
У › р т-1 р 


(рг” = анОВ Аг”, ты —@“ В" Аг”) —> ша › (28) 


т-1 
которая равносильна задаче 
(Рг”, =") -т (РВ Аа”, 2") (Р”,В Аа”) + 
+172 (РВ А2",В 142") —> п. 
Минимум последнего функционала достигается при значении итерационного параметра 
_ (РВ Аг", 2") +(Рг",ВТ Аг”) 
2(РВТ Аз", ВА") | 





(29) 


т-+1 


Дальнейшее зависит от выбора энергетического пространства Н р. В частности, итерационные методы вари- 
= 
ационного типа (27), (29), для которых р = А"В 'А ‚ приводят к следующим равенствам: 


— (ААВТАВ Аз", =") +(А"В Аг", В 'Аг”) _ 
2(А’ВТАВ Аз", В Аг”) 





т-1 


_ (АВТ Аз", В Аг") + (АВ А", В Аг”) _ 
2(В'АВ'Аз", АВ\Аг") 





_ (АВ 'Аг", В 'Аг")+(АВ'Аг", В Аг") _ (АВ'Аг", В Аг") 








2(В'АВ'Аз", АВ 'Аг”") (ВТАВ Аз", АВА"). 
Следовательно, итерационные параметры определяются ИЗ соотношений: 
(А»”, „”) 
Чи = 1 3 Ву" = Ах" — т т=0,1,... (30) 
(в` Ам”, Аи”) 


3. Исследование сходимости метода минимальных поправок. Если, как это принято, и" = В'Аз”" — век- 
тор поправки, а векторы погрешностей определены равенствами (26), то из уравнений (23), (24) получаем: 
Ва" = В" 1,142”. (31) 
Представим последнее соотношение в виде: 
И" он Ам". 
Выполним замену у” = В", С = В АВ!” и из последнего равенства получим: 
Вт В а рвут, 


что равносильно равенству 





























У у" аииСУ" = (Ев С)у”. (32) 
Используем оценку |” | 
т-1 т т 
У = |(Е-а„С)у 5 (ее -а„ыСо)+ (1-9) Евы, )}» | 
Выполняя очевидные преобразования, из последнего равенства получим С =Су+ С, ‚ где Су = и ‚= С" з 


Имеем очевидную оценку: 


ут 














ви [Е-аеась + (1-0, )Е-оиыС, м". (33) 


т-1 


Проведем оценку первого слагаемого в правой части неравенства (33): 
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Е— и С ут 
9, 


в [= ый = влад)” в = вил — 


т-1 т-+1 


2 


2 
| _ ин сут 


т+1 


























2 (34) 
= (Вт, ВУ" ) —2 В (дм”, и’ ) + [== (Ам", ВТ Ари" | 


т-1 


тот 
0-1 (Ах вы | 
Примем во внимание, что минимум дроби 0. +1 / 9. т достигается при = 


6 (В "Ам", Ади") 





Тогда оценку (34) запишем: 


2 
| — о 2 ут 
т+1 


ях 


В! Адм" , Арм”" 


быки) 
(В "Ал", Ам”) 














— (вт в" 2 (дм, )+ 








и] 
(В 'Ам", Ам") 








2 
1 (дм, в7"Ави” ) > (Ви, ВИЗ" ) 








ь 2 
СЕ °з 
т у") 


(с»", Сум" ) (С»" ‚Соу” ) ( : 


Таким образом, из (35) получаем равенство — оценку для первого слагаемого, стоящего в правой части нера- 


Е-— 9-1 а" ут 
Ви 


Из соотношения (33) для второго слагаемого с учетом равенства (см", у") =0 получаем: 


венства (33): 






о") 
(си", Сом" (> 














| в 


2. 
тот 
‚У 














(-е,)Е-а„ы@)у" | = (1-9.)>”, (1-6, =)у”)- 


т т т т 
-(в„„.Сь ‚(1-61 += (6 +1С Чи СУ )= 








2 0+1 
= (1 Фа} (Ут у" +02 Сл", 
ь ) ” В 
С учетом полученных неравенств (36), (37) из соотношения (33) имеем: 
(еше, 
0,1 
2 
(с»" , ут ) 


РЕ 


[0 [0 
Ът-1 бу", Ът-1 с") = 
0 


(1 —0,„.1 | (о ь" ) +0 0 


0 т 
ее (37) 
О 


у 


<60 





























Кое, и)Е-и,ыс)- 


т-1 


















т-1 т-+1 
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[02 [0 
т-1 Су у". т-1 С "| 
2 Ве Ва | "| 
——«„ КИ ||. 
тот 
(› ыы ) 


Введем обозначения: 





2 Чт т Я тт т 
(сь»",»") Е Су , 7. Су 


5 ^. 1 т-+1 т-1 


т- › Утшы= 
(сиу" сиу" ("м") ("у") 











(38) 


Определим итерационный параметр следующим образом: 
(ло” : о" в 
(В'Аю", Ак”) 





(39) 


0+1 


т-1 


Оценка для |У 














с учетом соотношений (38) и (39) примет вид: 





|" 








= (6. я /- 20. + 0 (1 +7) у) "|. 


г 1 Пин: 
1+ Ут 


| 


Положим, 0 . Тогда из последнего равенства получим оценку: 


т-а = 





т-1 


2 
< 1 — Пи- . Е Ут — 2+ 2 + (1 — Пи- ) 
= 1+ и т-1 


























1+ У т-=1 





-2+ 2+1 Бй 5 2+1 Е ЦА 
1+ Ут-ы 


1+ Ут 





1 Пин 5 
1+у 





т-1 














т-1 


Таким образом, последнее неравенство запишем в виде 





Ут ин + Пи 
1+ Ут 


ут < 1- Пи $ 


1 т-1 + 
+ Ут 














|". (0) 


Чтобы выбрать оптимальное значение параметра П„., вычислим производную от выражения, входящего в 


правую часть неравенства (40), и приравняем ее нулю: 


5 











: т-1 + Пи =0. (41) 
ы Ут (1 + У т-т ) [а + и ) 





т-1 


Заметим, что 


>0. 





сы Пи 


Тогда точка минимума выражения, входящего в неравенство (40), достигается при оптимальном значении па- 


раметра пи: 








2 
= ани. | (42) 
(1 У — бт ) 
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т-1 


у 


С учетом полученного оптимального значения Пи. , подставляя (42) в (40), получим оценку для | 


























2 2 
1— Эт Ут ЗУ тн 
1 2 Ут + 2 
тч | < Ут — ть тн т|| _ 
У < я бт + Я У 
+ У т-+1 + У т-+ 
Ут 
Зи — бт- 1 2 
№ Уж — би + Ут 
2 
1+ У т-ы 1+ Ут- — т-ы 





2 
бт + ти (1 Ут — бин ) 


1+ У т-н 





т 
у 














Константа р, определяющая скорость сходимости итерационного метода, оценивается сверху следующим 


К (| +7, — и 
= МЕРЫ. (43) 


1+ У т--1 


образом: 


Скорость сходимости метода минимальных поправок зависит от: 











Пусть Яма (С ) ‚ Апах (Со) — соответственно минимальное и максимальное собственные числа оператора Су, 


ау=А их (Со )/А ты (Со) — число обусловленности оператора Со. 


2 
Применим неравенство ху < (ах+ У/ а) /4 ‚ справедливое для а=0. 








Тогда 
ги <! ый ] < 
(е(лв Азии") + (Вит а) 
<1 - = 
(алиах (Со) +1 ани (Со)) 





Положим, а = а (С, )/ М (Со) . Тогда выражение, стоящее в правой части соотношения (43), меньше 


единицы при 5„.|<1. Для 5.1 И ум. имеют место оценки: 


т-+ т- 








= А (44) 
2 тах › 
(у +у?] У+1 


Введем обозначение 


—1 т т 
Киа = т ды”) ; Ува = Юнна (1-2) : (45) 
(В'Ам", Ави") 
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Выполним минимизацию оценки сверху постоянной р. Вычислим производную по 5„.| от выражения в пра- 


вой части (43) и, учитывая (45), получим: 


бт + (1 и ЗИ ) Аи (Кин + 1) 
1+ (1-52) 





Ят-1 


_ (1 + Ат )-2 Ут- Ат (а ЕО 1) + а 











- 20. (46) 
в 
(1 + Аи (- Ят-т )) 
Из неравенств (43) и (46) получаем оценку: 
2 
тах + и (1 + Ут — тах ) 
р< ь (47) 


1+ У тн 


2 
Выполняя тождественные преобразования — умножая на (о АУ (| ах ] / (1+7и-1) чис- 


литель и знаменатель выражения в правой части неравенства (47), получим: 


м + тины (1 + Ут — бах У (ин ть тн (1 + Ут-ы — о } + У т- ) 
р< — 
бтах Ут-Н + Ат (1 Уи — в ) 


(= (1 + У т-н ) + У тах а + 1) У т- (1 + Ут — гы } по У) 


5 тах Ут + 7» (1 о и ) 





В результате из неравенства 





У т-+ бя + ны (| Ут — т ) 




















Р < мах , 
2 
тах Ут + на (1 + Ут-н1 — бтах ) 
в силу (44) получим: 
Ут-ы (1 т Ут — г Е Ут-ы Ут-ы (| в Ут — га, + Ут 
р<|у 1 У +1 
52. 2 
У т-н (1 Уж — „)- Ут Ут (1 + Ут-1 — 5тах ) Ут 
(48) 
Введем обозначение 
е Ут-ы (1 + Ут- — Си + Ут-ы 
У =У , (49) 
52. 
Ут (1 в Ут — „)- Ут 
Тогда с учетом оценки (48) и введенного обозначения получаем 
= 
в (50) 
у +1 


Равенство (49) с учетом (45) примет вид: 


& К (1- ое )(1 А (1- А )- 52, м на (1 = ед ) 
У =У = 


А-а (:- га )(1+ Аи (| в )-52. я Аи (' ян) 


2 
(ЕК) +% (а) +6) ИИ 
он ыы) + т+1 _ н( н) - = 1+ Аи + . 


У 
Кин (1 + НЫ ) — ет К 
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Таким образом, в постоянную р, определяющую скорость сходимости метода минимальных поправок в соот- 
ветствии с неравенством (50), входит величина 


А И (1) 


где у — число обусловленности матрицы Су ‚а Ат характеризует отношение нормы кососимметрической части 


оператора А к симметрической: 


а [к 
т- [Аь” я и хо [Соу| 


На этом построение и оптимизация МАММП завершены. 


| 











Приведем результаты численных экспериментов, используя в качестве тестовой модельную задачу, приведен- 
нуюв [1]. В качестве характерного значения коэффициента сжимаемости возьмем В; = 10-3 /МПа . Остальные данные 
— такие же, как для тестовой задачи в статье [1]. В таблице приведены результаты численных экспериментов — ко- 
личества итераций, необходимые для достижения заданной точности & = 10° тремя методами: Зейделя (МЗ), верхней 


релаксации с шахматной упорядоченностью узлов (МВРШУ) и модифицированного адаптивного метода минималь- 
ных поправок (МАММП). 























Таблица 1 
Количества итераций, требуемые для решения модельной задачи употребительными методами 
Номер временного слоя Метод решения сеточных уравнений для функции давления 
МЗ МВРШУ МАММП 
1 32535 7342 2752 
2 20347 3453 1014 
3 13138 2387 696 
4 12190 2123 688 
э 12132 2135 692 




















Заключение. В статье описано построение и исследование модифицированного адаптивного метода мини- 
мальных поправок — МАММП, предназначенного для численного моделирования пространственно-трехмерных про- 
цессов фильтрации двухфазной сжимаемой жидкости в пористых средах. Такие процессы являются типичными в за- 
дачах моделирования нефтяных месторождений — в случае, когда нефтеводоносный пласт находится под воздействи- 
ем давления, изменяющегося значительным образом, и необходимо учитывать сжимаемость среды. Для данного клас- 
са сеточных задач с несамосопряженным оператором построен МАММП. Проведена его оптимизация и доказана схо- 
димость метода, который обладает асимптотической скоростью сходимости, характерной для «классического» ПТМ, 
не использующего технику чебышевского ускорения и применяемого для задач с самосопряженным оператором. Ре- 
зультаты численных экспериментов продемонстрировали эффективность метода. Число итераций уменьшилось в 3-20 
раз по сравнению с другими часто применяемыми подходами к решению сеточных задач фильтрации (методом Зейде- 
ля, методом верхней релаксации и другими). 


Библиографический список 


1. Сухинов, А. И. Усовершенствованный попеременно-треугольный метод численного решения простран- 
ственно-трехмерных задач фильтрации двухфазной несжимаемой жидкости / А. И. Сухинов, Л. А. Григорян, 
А. А. Сухинов // Вестник Дон. гос. техн. ун-та. — 2016. — Т. 16, № 1. — С. 5-18. 

2. Коновалов, А. Н. Задачи фильтрации многофазной несжимаемой жидкости / А. Н. Коновалов. — Новоси- 
бирск : Наука, 1988. — 166 с. 

3. Вабищевич, П. Н. Явно-неявные вычислительные алгоритмы для задач многофазной фильтрации / 
П. Н. Вабищевич // Математеческое моделирование. — 2010. — Т. 22, № 4. — 118-128. 

4. МафетаНса] то4деПп» оЁ Но\’з ш рогоиз шеФа / У. В. изба [её а] /УУЗЕА$ Тгапзасйоп$ оп ИА тесВап- 
1с$. — 2014. — Ус. 9. —Р. 130-166 





Сухинов А. И. и др. Математическое моделирование фильтрации двухфазной сжимаемой жидкости 


5. ЗаякКШтоу, А. 1. Адарнуе шо е4 аКегпаНп® шлапоШаг Цегануе тефо4 Гог зо]уше 14 едиаНоп$ \И а поп- 
зе{-аЧое орегаог/ А. 1. ЗаКЫтох, А. Е. СызуаКоу // Мафетайса! Мо4е|5 апа Сотршег ЗпишаНоп$. — 2012. — 
Уо/. 4, 133. 4. — Р. 398—409. 

6. ЗаяКШтох, А. 1. шсгеазше еЁЯслепсу оЁ а№егпайп® шапоаг тефо4 Ъазе4 оп пиргоуе зресша] езИтлавез / 
А. 1. ЗаКыпох, А. У. 5ЫзБепуа // МафетаНса! Мо4е[5 апа Сотрщег ЗииаНоп$. — 2013. — У\о1. 5, 155. 3. —Р. 257—265. 

7. Сухинов, А. И. Модификация метода минимальных поправок для решения сеточных уравнений с несамо- 
сопряженным оператором / А. И. Сухинов, А. Е. Чистяков, А. В. Шишеня // Известия ЮФУ. Технические науки. — 
2013. — № 4. — С. 194—202. 





ВеГегепсе$ 


1. баКЫтоу, А.Т., Сиеогуап, Г.А., ЗаКЫтпоу, А.А. ОзоуегзВепуоуаппуу рорегетеппо-неиеоГпуу теюа 
сыЗеппого гезвешуа ргозбгап$уеппо-шекытегпукВ тадась Н“гази ауиКОатпоу пез7б1таетоу 75 Кози. [Адуапсе4 а1- 
{егпае лапошаг те#о4 о питепса] зоайоп о зраНа| гее-Читепз1опа! ргоеглз оЁ смо-рБазе Яйганоп оЁ шсотргезЫе 
Ниа.] УезциК оЁ О$ТО, 2016, уо1. 16, по. 1, рр. 5—18 (т Киззлап). 

2. Копоуа]оу, А.М. Хадасы ЯГгави тшпозоатпоу пезтВитаетоу 751АКози. [Ргоеглз оЁ по!арБазе Нганоп оЁ 
шсотргезЫе Яла.] Моуозизк: Мака, 1988, 166 р. (шт Виззап). 

3. УабзЬсВеу1сь, Р.М. Уаупо-пеуаупуе уусб15ШеГпуе а!гопту Фуа хадасВ тпозоа7поу Нави. [ЕхрПси- 
парНсй питепса| а]готИт$ ог шыа-рБазе ВИтайоп ргоЫетз.] МафетаНса1 Мо4е]5 ап Сотрщег Зпищайопз, 2010, 
уо1. 22, по. 4, рр. 118—128 (шт Визчап). 

4. РизЬа, У.В., её а1. Мафетайса| тодеНпе оЁ Ио\з ш рогойз те фа. \УЗЕА$ ТгапзасНопз$ оп ИА тесВашсз, 
2014, ус. 9, рр. 130-166. 

5. баякытоу, А.Г., Сызуакоу, А.Е. Адарнуе тод1еа аПегтпайпе шапочаг иегануе тефо4 Гог зо]ушз 214 едиа- 
00105 УЛ а поп-зе{-ао1те орегаюг. Мафетайса! Мо4е!5 апа Сотршег ЗпишаНопз, 2012, уо/. 4, 133. 4, рр. 398—409. 

6. ЗякКытох, А.Т., 5515Вепуа, А.У. псгеазте еЁйсепсу оЁ а№егпайпе шлапоу]аг тейо4 Базе оп ппргоуе4 5ресёга1 
езитаез. Мафетайса! Мо4е!$ ап Сотрщег Зиаайопз, 2013, уо1. 5, 155. 3, рр. 257-265. 

7. Закытоу, А.Т., СызбуаКоу, А.Е., ЗшзВепуа, А.У. Моа!кКазтуа теюда пипипаГпукВ роргауоК Чуа гезвешуа 
зеосбпуКВ игаупешу $ пезатозоргуатВеппут орегабюгот. [Мо@1ЙсаНоп оЁ плшита]| гез14иа]$ Иегануе тефо4 Рог зоУтз 
эт14 едиаНоп$ УИ поп-зеа ош орегагогз.] [7уезнуа ЗЕе40. Епэзтеенпе Зс1епсез, 2013, по. 4, рр. 194—202 (ш ВКаззап). 


Поступила в редакцию 19.06.2016 
Сдана в редакцию 20.06.2016 
Запланирована в номер 07.07.2016 


